Fate of the first traversible wormhole: black-hole collapse or inflationary expansion 
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We study numerically the stability of Morris & Thorne’s first traversible wormhole, shown previ- 
ously by Ellis to be a solution for a massless ghost Klein-Gordon field. Our code uses a dual-null 
formulation for spherically symmetric space-time integration, and the numerical range covers both 
universes connected by the wormhole. We observe that the wormhole is unstable against Gaussian 
pulses in either exotic or normal massless Klein-Gordon fields. The wormhole throat suffers a bifur- 
cation of horizons and either explodes to form an inflationary universe or collapses to a black hole, if 
the total input energy is respectively negative or positive. As the perturbations become small in to- 
tal energy, there is evidence for critical solutions with a certain black-hole mass or Hubble constant. 
The collapse time is related to the initial energy with an apparently universal critical exponent. For 
normal matter, such as a traveller traversing the wormhole, collapse to a black hole always results. 
However, carefully balanced additional ghost radiation can maintain the wormhole for a limited 
time. The black-hole formation from a traversible wormhole confirms the recently proposed duality 
between them. The inflationary case provides a mechanism for inflating, to macroscopic size, a 
Planck-sized wormhole formed in space-time foam. 
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I. INTRODUCTION 


Wormholes are known as a kind of solution of the 
Einstein equations, and have become a popular research 
topic, raising theoretical possibilities of rapid interstel- 
lar travel, time machines and warp drives. These topics 
sound like science fiction, but after the influential study 
of traversible wormholes by Morris & Thorne [1], it be- 
came widely accepted as a scientific topic [2]. The only 
physically non-standard feature is that one has to assume 
negative-energy matter to construct such a wormhole. 
However, such exotic matter occurs extensively in quan- 
tum field theory and in alternative gravitational theories 
such as scalar-tensor theories. 

Until recently, wormholes have been studied mainly 
as static or cut-and-paste models, or without an inde- 
pendently defined exotic matter model. Not long ago, 
one of the authors [3] proposed a unified theory of black 
holes and traversible wormholes, arguing that the two are 
dynamically interconvertible, and that traversible worm- 
holes are understandable as black holes under negative 
energy density. This opens a new viewpoint on the dy- 
namical nature of both black holes and wormholes, in- 
cluding Hawking radiation. This synthesis has been ex- 
amined using a low-dimensional model [4], where the the- 
ory is affirmed. 

The purpose of this article is to investigate wormhole 
dynamics in four-dimensional Einstein gravity, using nu- 
merical simulations. Our starting point is a static worm- 


hole which is perhaps best known as Morris & Thorne’s 
opening example (their box 2). This metric is actually a 
solution for a massless Klein-Gordon field whose gravita- 
tional coupling takes the opposite sign to normal, as was 
shown earlier by several authors [5], the earliest of which 
appears to be Ellis [6], who called it a drainhole. To 
our knowledge, this is the earliest solution which would 
nowadays be called a traversible wormhole. 


We study dynamical perturbations of this static worm- 
hole, using the spherically symmetric Einstein system 
with the above exotic matter model, the massless ghost 
Klein-Gordon field. We developed a numerical code 
based on a dual-null coordinate system, in order to follow 
the horizon dynamics and radiation propagation clearly. 
Our main experiment is to add or subtract Gaussian 
pulses in the ghost field, i.e. respectively with negative 
or positive energy. We also consider Gaussian pulses in 
a normal Klein-Gordon field, to see the effect on the 
wormhole of normal matter, like a human being travers- 
ing the wormhole. We discover how the initially static 
wormhole will change its structure due to these dynamic 
perturbations. Although our model is specific to spheri- 
cally symmetric space-times, we believe that it illustrates 
the dynamical nature of traversible wormholes. To our 
knowledge, this is the first numerical study of wormhole 
dynamics. 


We describe our model and numerical method in $I 
and present numerical results in §III. §IV concludes. 


Il. MODEL AND NUMERICAL METHOD 
A. Field equations 


The field equations for a massless conventional and 
ghost Klein-Gordon field, yY and ¢ respectively, in Ein- 
stein gravity are, in standard notation, 


R= 2Vyp@ VY — 2768 Vd, 
O¢ = 0, Oy =0 


(2.1) 
(2.2) 


where the coupling constants have been fixed by choice 
of units for y and @. We adopt a dual-null formulation 
developed by one of the authors [7]. The spherically sym- 
metric line-element may be written in dual-null form 


ds? = r2dS? — 2e~fdatdxr- (2.3) 


where dS? refers to the unit sphere. Then the configu- 
ration fields are (r, f,¢,) as functions of (x+, £7). The 
Einstein equations in spherical symmetry are expressed 
in [8, 9], leading to the field equations 


0.047 + (Os f)(Ozr) 

= —r(O,~)? + r(O4¢)?, 
r?04,0_f +2(O,r)(O_r) + ef 

= 2r? (ð Y) (3- Y) — 2r? (34H) (3-4), 
r0,0_r + (0,r)(0_r) +e-F/2 = 0, 
r0,0_¢ + (4r)(3-¢) + (3-r)(3+¢) = 0, 
r040_wv + (04r)(O_W) + (O_r) (04%) = 0. 


(2.4) 


2.5) 
2.6) 
2.7) 
2.8) 


~~ oa > 


To obtain a system accurate near S*, we can use the 
conformal factor Q = 1/r instead of r and, as first-order 
variables, the conformally rescaled momenta 


04 = 24r = —-2077040, (2.9) 
ve = Oxf, (2.10) 
p+ = rô = Ql, (2.11) 
Te = rep =O" dz. (2.12) 


Here + has the meaning of expansion, though it is 
rescaled from the usual definition, 04 = 2r~!Ozr. Then 
(P+, T+, V+, V+) are finite and generally non-zero at SF. 
The remaining first-order equations follow from the above 
field equations and the spherically symmetric identity 


0- = O_O4: 


0404 = -v4 —20n2 + 2007, 2.13) 
ðv = —A(9,0_/2+e7%), (2.14) 
zvy = —07(0,0_/2+e7 

—2nr47_ +2949_), (2.15) 
tp = —WWz04/2, (2.16) 
dite = —NW-z74/2. (2.17) 


These equations plus the inverted momentum definitions 


049 = —0794/2, (2.18) 


anf = v4, (2.19) 
ðo = Ops, (2.20) 
Os = Arg (2.21) 


constitute the first-order dual-null form, suitable for nu- 
merical coding. 


B. Numerical method 


We prepare our numerical integration range as drawn 
in Fig.1. The grid will cover both universes connected by 
the wormhole throat xt = x~. The grid is assumed to 
be equally spaced in Az = A. 

The basic idea of numerical integration is as follows. 
We give initial data on a surface S and the two null hy- 
persurfaces X4 generated from it. Generally the initial 
data have to be given as 


(Q, f v+, Q, p) 


(Vi, 04,74) on Ne: x 


We then evolve the data u = (0,04, f, v+, Q, Y, p+, T+) 
on a constant-a~ slice to the next. 

Due to the dual-null decomposition, the causal re- 
gion of a grid is clear, and there are in-built accuracy 
checks: the integrability conditions or consistency con- 
ditions 0_O,u = 0,0_u. In order to update a point N 
(north), we have two routes from the points E (east) and 
W (west). The set of equations, (2.13)-(2.21), gives us 
x*-direction (W to N) and «~-direction (E to N) inte- 
grations together with the consistency conditions. 

All the equations can be written as 


Unew = Uold + f(urms)A, (2.24) 


and the input variables, urpas, in the right-hand-side 
term depend on the integration schemes. In more detail, 
we made the steps: 


1. Integrate (Q, 94, 0_,v_, f,¢,¥, @_,7_) using the 
-equations (update N from W): use urns = uw. 


2. Integrate (0,04, 0_,74,f,¢,¥, 04,74) using the 
O_-equations (update N from E): use urns = ug. 


3. Update N from W again: use urns = (uwtun)/2. 
4. Update N from E again: use uggs = (ug +un)/2. 


5. Check the consistency of (0,01, 0_, f, p) at point 
N. If the convergence is not satisfactory, repeat 
back to the step 3. We set the tolerance as 
max \un—+N = uws] < 107° for (Q,04,0_, Í, Q). 


As a virtue of the dual-null scheme, we can follow the 
wormhole throat or black-hole horizons easily. They are 
both trapping horizons, hypersurfaces where J+ = 0 or 
V = 0 [9, 10]. Another benefit is the singular point exci- 
sion technique. As we described, the causal region of each 


grid point in the dual-null scheme is apparent. When a 
grid point is inside a black-hole horizon and near to the 
singularity, we can exclude that point and grid points in 
its future null cone from further numerical computation. 
Actually we excised grid points if one of (V4 |, |Y_|, |e~“]) 
exceeds ~ 25. We found that this is useful for investigat- 
ing inside a black-hole horizon, though there is naturally 
a limit on how closely a singularity can be approached. 


C. Static wormhole 


The first Morris-Thorne wormhole metric [1] was given 
as 

ds? = (a? + 1°)dS? + dl? — dt? (2.25) 

where a denotes the throat radius of the wormhole. Since 

it is an overall scale, we set a = 1 for numerical purposes, 

but retain it in the text. Transforming the proper radial 


length l and the static proper time t to dual-null coordi- 
nates 


at = (t+l)/v2 (2.26) 


the analytic expressions for the Ellis wormhole solution 
[6] are 


@ = tan™!(l/a), (2.27) 
Q = 1/Va2 +P, (2.28) 
f = 0, (2.29) 
ps = ta/V2Ve2 +22, (2.30) 
84 = +V2i/Vo2 +P, (2.31) 
vs = 0 (2.32) 


with the conventional field ~ vanishing. The wormhole 
is symmetric in interchange of the two universes, | + —1, 
though this will not be so for the perturbations. The 
energy density is 


a2 


(a 
i 8r (a? 4+ 17)? 


(2.33) 


with units G = 1. The negative energy density is a max- 
imum at the throat l = 0, tends to zero at infinity and 
is smaller for larger wormholes. Note that the wormhole 
throat is a double trapping horizon, 0+ = J_ = 0. The 
behaviour of these generally different trapping horizons 
will be a key indicator of the evolution of the perturbed 
wormhole. 
The dual-null initial data for the static wormhole is 


(0, Q, f, 0+) = (0, 1/a, 0,0) on S$ 


(94,V4) = (ta/\/2a? + (a*)?,0) on H+(2.35) 


D. Gravitational mass-energy 


We find that the dynamical structure will be charac- 
terized by the total mass or energy of the system, the 
Bondi energy. Localizing, the local gravitational mass- 
energy is given by the Misner-Sharp energy E, while the 
(localized Bondi) conformal flux vector components y* 
were defined in [8, 9]: 


E = (1/2)r[1 — g`! (dr, dr)] = (1/2)r + efr (O,r)(0_r) 


1 1 
= — [1 + że 5 2: 
5a + 5e 040_] (2.36) 
gt = r°Tttéar = ref T +ô4r 
= e” (nÈ — p )04/8r. (2.37) 


They are related by the energy propagation equations or 
unified first law [9] 


OLE = 4rps = -iei — pox, (2.38) 


so that one may integrate to 


(at ,2~) 
RG er) = stan | (p,dx++p_dx~), (2.39) 
2 (0,0) 


where the integral is independent of path, by conserva- 
tion of energy [8, 9]. Here the initial condition E|s = 
a/2 corresponds to the static wormhole. The defini- 
tions are local, but lim,+_,,., E is the Bondi energy and 
lim,+ 35. y— the Bondi flux for the right-hand universe. 
For the static wormhole, the energy E = a?/2Va? + l? 
is everywhere positive, maximal at the throat and zero 
at infinity, 1 + -too, i.e. the Bondi energy is zero. Gen- 
erally, the Bondi energy-loss property, that it should be 
non-increasing for matter satisfying the null energy con- 
dition [8], is reversed for the ghost field. 


Ill. DYNAMICAL BEHAVIOUR OF 
TRAVERSIBLE WORMHOLES 


A. Evolution of static wormhole 


We start with the numerical evolution of the static 
wormhole. We find that numerical truncation error can 
quite easily destroy the stability. In Fig.2, we show our 
tracking of the location of the trapping horizon ð- = 0 
for static wormhole initial data, with different numbers 
of grid points: 801, 1601, ... 9601, for xt = [0,20]. The 
location of the throat, where Ji = V- = 0, is supposed 
to be stationary at xt = x~. The figure shows that 
at sufficient resolution, we can maintain a static worm- 
hole for a given time. Vice versa, this shows that in our 
dual-null scheme and the current applied numerical inte- 
gration scheme, the lack of resolution inevitably causes 
an unstable evolution eventually. We hereafter discuss 
the dynamical structure of perturbations of the wormhole 


wormhole throat 


FIG. 1: Numerical grid structure. Initial data are given on 
null hypersurfaces X+ (a+ = 0, x* > 0) and their intersection 
S. 
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FIG. 2: Convergence behaviour of the code for exact static 
wormhole initial data. The location of the trapping horizon 
V = 0 is plotted for several resolutions labelled by the num- 
ber of grid points for «+ = [0,20]. We see that numerical 


truncation error eventually destroys the static configuration. 


only in the reliable evolution range, with the resolution 
kept good enough (9601) to preserve the static case. 

In Fig.3, we demonstrate the numerical evolution of the 
static wormhole, giving the expansion 0, and the local 
gravitational mass-energy E as functions of (£t, x7). 


B. Gaussian perturbations in ghost field 


We first put perturbations of the static wormhole in 
the form of Gaussian pulses, input from the right-hand 
universe l > 0. The initial data on X+ is 


p+ = a/V2/ a? +1? + ca exp(—cx(l— c)”), (3-1) 


expansion plus 
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FIG. 3: Static wormhole configuration obtained with the 
highest resolution calculation: (a) expansion 3+ and (b) lo- 
cal gravitational mass-energy E are plotted as functions of 
(xt,x7). Note that the energy is positive and tends to zero 
at infinity. 


with all other initial data as for the static wormhole. Here 
Ca; Cb; Ce are parameters, with runs being performed for 
increasingly small values of the amplitude cg, cp = 3,6,9 
and ce = 3. That is, the pulses will hit the wormhole 
throat at xt = x = 3. For Ca positive or negative, 
this corresponds respectively to enhancing or reducing 
the supporting ghost field. Note that this pulse-type per- 
turbation starts only from the right-hand universe. This 
situation is available because our numerical region covers 
both sides of the wormhole, though we assume that the 
space-time is locally spherically symmetric. 

Fig.4 shows the horizon locations for both J, = 0 and 
v_ = 0, for cy = 3. We plotted three cases, cg = £0.1 
and —0.01, and it can be seen that in each case, the 
double trapping horizon comprising the initial wormhole 
throat bifurcates when the pulse hits it. The two horizons 
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FIG. 4: Horizon locations, V+ = 0, for perturbed wormhole. 
Fig.(a) is the case we supplement the ghost field, ca = 0.1, and 
(b1) and (b2) are where we reduce the field, ca = —0.1 and 
—0.01. Dashed lines and solid lines are 0+ = 0 and v- = 0 
respectively. In all cases, the pulse hits the wormhole throat 
at (x", x7) = (3,3). A 45° counterclockwise rotation of the 
figure corresponds to a partial Penrose diagram. 


then accelerate away from each other, rapidly approach- 
ing the speed of light and appearing to be asymptotically 
null. Due to the grid excision technique for strong cur- 
vature regions, the final horizon coordinate xp, where 
the right-hand trapping horizon becomes null, is often 
excised from the numerical computation. In such a case, 
we determine zy by extrapolating the horizon trajectory 
with a function xt = bo +bia~ +b2(a~)? +b3/(x~ — £p), 
where b; and xh are unknown constants. A partial Pen- 
rose diagram of the evolved space-time is given in Fig.5, 
which is like the prediction in [3]. However, the two 
horizons move in opposite directions depending on the 
sign of Ca. In the following, one should bear in mind 
that J, > 0,v_ < 0 initially in the right-hand universe 
x* > x, with signs reversing for the left-hand universe. 

When we supplement the ghost field, Ca = 0.1, we see 
that the v+ = 0 horizons move respectively in the x? 
directions, meaning that J > 0 in the region between 
them. This defines past trapped surfaces, which occur 
in expanding universes, such as the de Sitter solution. 
In Fig.6 we plot the areal radius r of the “throat” [17] 
xt = x7, as a function of proper time 


r= | Ea t=(et+a-)/V2. (3.2) 


It can be seen to rapidly increase, actually exponentially: 
the lines in the figure can be expressed with a function 
r/a = 1+ bsexp(H(7 — b5)) where b; and H are con- 
stants, and we show the Hubble constant H in Table I. 
This exponential expansion, combined with the horizon 
structure, indicates that the wormhole has exploded to 
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(b2) pulse input with smaller positive energy 
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FIG. 4: (continued) 


an inflationary universe. The two trapping horizons be- 
come the cosmological horizons. This connection be- 
tween wormholes and inflation was unanticipated, but 
is consistent with the general theory [3]. The two uni- 
verses connected by the wormhole have essentially been 
combined into one universe. 

On the other hand, when we reduce the ghost field, 
Ca = —0.1, the V+ = 0 horizons move respectively in the 
a* directions, meaning that J, < 0 in the region be- 
tween them. This defines future trapped surfaces, which 
occur inside black holes, such as the Schwarzschild solu- 
tion. Fig.6 shows that the areal radius is heading towards 
zero in finite proper time, though the code cannot follow 
it into the apparent central singularity. This and the 
horizon structure indicates that the wormhole has col- 
lapsed to a black hole. The trapping horizons appear to 
become asymptotically null, where they would coincide 
with the event horizons of the black hole. This occurs 
around xp ~ 4.46 for the right-hand universe. 

If the amplitude of the pulse is reduced, the black hole 
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FIG. 5: Partial Penrose diagram of the evolved space-time. 
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FIG. 6: Areal radius r of the “throat” «* = x`, plotted as 


a function of proper time. Additional negative energy causes 
inflationary expansion, while reduced negative energy causes 
collapse to a black hole and central singularity. 


forms later, as shown in Fig.4(b2) for ca = —0.01. As a 
physical measure of the size of the perturbation, we take 
the initial Bondi energy 


Eo = E(%,0) (3.3) 


scaled by the initial throat radius a. In practice we take 
E(20,0) and subtract the corresponding value for the 
static wormhole. Some data are presented in Table I. We 
find that positive Æo will cause collapse to a black hole, 
while negative Eo will cause explosion to an inflationary 
universe. Also, the speed of the horizon bifurcation in- 
creases with the energy. One could also scale the total 
energy by the maximal initial energy E, but this occurs 
at the throat and is a/2. In either case, the perturbations 
are small, down to 1% in energy, yet the final structure is 
dramatically different. Thus we conclude that the static 
wormhole is unstable. 

Accepting that the static wormhole is stable to lin- 
ear perturbations [11, 12], we seem to have discovered a 
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(b1) pulse input with positive energy 
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(b2) pulse input with smaller positive energy 
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FIG. 7: Energy E(at,2~) as a function of a7, for at = 


12, 16,20. Here ca is (a) 0.05, (b1) —0.1 and (b2) —0.01. The 
energy for different «* coincides at the final horizon location 
£p, indicating that the horizon quickly attains constant mass 
M = E(oo, xp). This is the final mass of the black hole. The 
values are shown in Table I. 


non-linear instability, something of interest in itself. Cir- 
cumstantial evidence for a quadratic instability is that it 
is first signalled by non-zero values of v+, and the prop- 
agation equations (2.15) for v+ have quadratic terms on 
the right-hand side, which should cancel to zero in the 
static solution. This could be addressed by second-order 
perturbation theory, but we know of no such studies for 
wormbholes. 


In Fig.7, we plot the energy E(x*,x~) as a function 
of £7, for x* = 12,16, 20 for the three cases in Fig.4. In 


(a) 


FIG. 8: Relation between the initial perturbation and the 
final mass of the black hole. (a) The final horizon (V+ = 
0) coordinate xp — 3, since we fixed ce = 3, versus initial 
energy of the perturbation, Eo. We plotted the results of the 
runs of ca = 1071,---,10~* with cy = 3,6, and 9. They lie 
close to one line. (b) The final black hole mass M for the 
same examples. We see that M appears to reach a non-zero 
minimum for small perturbations. 


each case, the mass increases rapidly in x7 as the horizon 
V+ = 0 forms, but rapidly approaches a constant value 
in t at the horizon. This value 


M = E(%, xp) (3.4) 
is the final mass of the black hole. Obtaining this easily 
is another virtue of the dual null scheme. The values are 
shown in Table I. The graphs indicate that an observer 
at infinity will see a burst of radiation as the wormhole 
collapses or explodes. For collapse, a certain fraction of 
the field energy radiates away, the rest being captured 
by the black hole, constituting its mass. For explosion, 
the radiated energy continues to rise in an apparently 
self-supporting way as the universe inflates. 

Fig.8 gives a survey of various parameters of c, and cp, 
fixing ce = 3. Fig.8(a) shows how quickly the black-hole 
horizon, the outer trapping horizon J; = 0, develops 
after the pulse hits the throat. We see that the dots 
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FIG. 9: Evolution of a wormhole perturbed by a normal 
scalar field. Horizon locations: dashed lines and solid lines 
are V4 = 0 and v- = 0 respectively. 


are close to one line, indicating a logarithmic relation 
between the collapse time and the initial energy of the 
perturbation, 

en ec, oe (3.5) 
Fig.8(b) is the final black hole mass M as a function 
of the initial energy of the perturbation. Interestingly, 


as the perturbations become small, the black-hole mass 
seems to approach a non-zero constant, 


M ~ 0.30a. (3.6) 
In the inflationary case, the Hubble parameter H (or 
horizon mass 1/2H) also seems to have a limiting value 
for small perturbations, as seen in Table I: 

H ~ 1.1/a. (3.7) 
Since these masses are the same order as the wormhole 
scale a, we guess that the static wormhole is something 
like an instanton with two different attractors nearby. 
This scenario requires more detailed study with different 
models, and we plan to report it subsequently. 


C. Gaussian pulses in conventional field 


Similarly, we next consider adding a small amount 
of conventional scalar field to the static wormhole so- 
lution. This is a simplified model of a situation that a 
positive-energy creature jumps into a traversible worm- 
hole. Again we add Gaussian pulses to the initial data 
on 4 as 


T} = Ča exp(—@(I — é)?), (3.8) 


(a) 


FIG. 10: The same plots as in Fig.8 for conventional field 
pulses. (a) The final horizon (v+ = 0) coordinate xj, — 3 
versus initial energy of the perturbation, Ho. We plotted the 
results of the runs of ča = 0.5,---,107? with & = 3,6, and 
9. They lie close to one line, different to that of Fig.8(a), but 
with the same gradient. (b) The final black hole mass M for 
the same examples. 


with all other initial data as for the static wormhole. We 
set Cy = Če = 3, so that the pulse will again hit the 
wormhole throat at 2t = z7 = 3. 

We show the horizon structure in Fig.9 and initial total 
energy in Table II. The additional normal field slightly 
enhances the total energy, with either sign of amplitude, 
and the wormhole collapses to a black hole. The fun- 
damental dynamical behaviour is the same as for the 
previous ghost perturbations with positive energy. We 
see again a logarithmic relation between the time and 
the initial energy of the perturbation, e?#~? œ ies 
Fig.10(a). The critical exponent is consistent with that 
found for perturbations in the ghost field, indicating a 
new kind of critical behaviour reminiscent of that found 
by Choptuik [13] in black-hole collapse. In the limit of 
small perturbations, the black-hole mass also appears 
to approach the same non-zero constant, M ~ 0.30a, 
Fig.10(b). Although the current code has a limit on how 
small perturbations can be, since the black-hole collapse 
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FIG. 11: Temporary wormhole maintenance. After a normal 
scalar pulse representing a traveller, we beamed in an addi- 
tional ghost pulse to extend the life of the wormhole. Horizon 
locations J, = 0 are plotted for three cases: (A) no mainte- 
nance, which results in a black hole; (B) with a maintenance 
pulse which results in an inflationary expansion; (C) with a 
more finely tuned maintenance pulse, which keeps the static 
structure up to the end of the range. 


should occur within the reliable evolution range, the data 
suggest the existence of a critical solution with a certain 
black-hole mass. 

For these perturbations, the energy perturbation is as 
small as about 0.01%, a quite convincing indication of 
instability. From this and general principles [3], we ex- 
pect that any normal matter traversing the wormhole will 
cause its collapse to a black hole. 


D. Wormhole maintenance 


Supposing the normal field pulse represents an actual 
traveller of the traversible wormhole, then he, she or it 
may go through the wormhole and exit safely into the 
other universe if the speed is high enough, as can be seen 
from the Penrose diagram obtained from Fig.9. However 
the wormhole itself will collapse to a black hole, ending 
its usefulness for travellers. 

We here demonstrate a kind of temporary mainte- 
nance of the wormhole, by sending in an additional 
ghost pulse just after the passing of the traveller. 
In Fig.11 we show the results for pulse parameters 
(Ga, Čb, Če) = (0.1,6.0,2.0) for the normal field repre- 
senting the traveller, combined with a balancing pulse 
(Ca; Cb, Cc) = (0.02390, 6.0, 3.0), case B, or (Ca, Cb, Cc) = 
(0.02385, 6.0, 3.0), case C. If we do not send the balanc- 
ing pulse, the wormhole collapses to a black hole with 
horizons given by case A in the plot. The case B ends up 
with an inflationary expansion, while the case C keeps 


final structure 


+10~'|3.0]—5.07 107° [inflation 1.56 (Fig.4a) 
+107?|3.0|—4.45 107° | inflation 1.19 
+107?|3.0|—4.40 1074| inflation 1.13 
+10~*]3.0] —4.39 10-5 |inflation 1.12 


>|black hole 8.42 0.296 
black hole 6.98 0.310 
black hole 5.64 0.325 (Fig.4b2) 
?|black hole 4.46 0.420 (Fig.4b1) 
black hole 4.82 0.38 
black hole 4.90 0.35 


t/a M/a 


TABLE I: Initial Bondi energy Eo and the final structures for 
different ghost pulse inputs described by (3.1). The energy 
is evaluated at xt = 20, subtracting that of the exact static 
solution. The final structure was judged at x™ = 10 by the 
horizon structure. For inflation cases, the Hubble constant H 
was measured on the “throat”, Fig.5. For black-hole cases, 
we show the final horizon coordinate xp, where the 94 = 0 
trapping horizon becomes null, Fig.4(b), and the final mass 
M of the black hole, Fig.7(b). 


t/a M/a 
-6.86 x 107°|black hole 9.61 0.31 
-1.71 x 107°|black hole 7.50 0.32 
+6.86 x 107*|black hole 6.68 0.34 (Fig.9) 
-1.71 x 107'|black hole 4.82 0.46 
+6.86 x 107} black hole 3.70 0.90 


TABLE II: Initial Bondi energy Eo and the final structures 
for different conventional pulse inputs described by (3.8). The 
columns are the same as for Table II. 


the wormhole structure at least until <~ = 10. Since 
the final fate of the wormhole is either a black hole or 
an inflationary expansion, to keep the throat as it was 
requires a fine-tuning of the parameters, and may not be 
realistic. However, it shows how the wormhole life may 
be extended. This indicates that the wormhole might 
be maintained by continual adjustments to the radiation 
level, though it would be a never-ending project. 


IV. DISCUSSION 


We have numerically studied the dynamic stability of 
the apparently earliest [6] and most frequently rediscov- 
ered [3, 5] traversible wormhole, most famous as the first 
Morris-Thorne wormhole [1]. We observe that the worm- 
hole is unstable against Gaussian pulses in both exotic 
and normal massless Klein-Gordon fields. The wormhole 
throat suffers a bifurcation of horizons and will either ex- 
plode to an inflationary universe or collapse to a black 


hole, if the total energy is respectively negative or posi- 
tive. For the normal matter, the total energy is necessar- 
ily positive and a black hole always results. This black- 
hole formation from a traversible wormhole confirms the 
unified theory of both [3]. 

The inflationary expansion provides a mechanism for 
enlarging a quantum wormhole to macroscopic size. 
Wheeler introduced wormholes to gravitational physics 
as part of his still-popular idea of space-time foam [14], 
which envisages transient wormholes at the Planck scale, 
due to quantum fluctuations in space-time topology. 
Morris & Thorne imagined pulling a wormhole out of 
this foam and enlarging it somehow. Actual mechanisms 
are scarce: Roman [15] suggested inflation, though with- 
out an independently defined exotic matter model, and 
Redmount & Suen [16] found a quantum-mechanical in- 
stability of certain cut-and-paste wormholes in a space- 
time foam model, also leading to unbounded expansion. 
We have discovered that for one of the simplest classi- 
cal exotic matter fields, wormholes can naturally inflate. 
The problem is only to stop the inflation proceeding in- 
definitely. 

The results also provide an entertaining answer to the 
question of what happens if someone attempts to traverse 
the wormhole. At best, with sufficient alacrity, our hero 
could exit into the other universe, only to look behind and 
see that the passage has caused the wormhole to collapse 
to a black hole, thereby sealing off the causal connection 
to the home universe. Perhaps there was also an unlucky 
companion who did not make it through, suffering the 
well-known grisly fate of being trapped in a black hole. 
The survivor would be left to ruminate on the ironic fate 
of those who would tamper with the fabric of space-time: 
stranded forever in the twilight zone. As we showed, 
such a tragedy could be avoided by a carefully calculated 
burst of ghost radiation, just before or after the traveller 
passes, which maintains the static wormhole for a certain 
time. The traveller then has a limited time to explore the 
other universe before the wormhole collapses. 

For our colleagues who are understandably sceptical 
about wormholes and exotic matter, we hope that our 
numerical study has at least demonstrated that they can 
be studied and understood by the same local, dynamical 
methods as used for black holes with normal matter. In- 
deed, the results lead straight back to mainstream physics 
of black holes and inflation, with new discoveries about 
both, particularly the apparent critical behaviour. For 
a general audience, perhaps the message is that (albeit 
theoretical) science can be stranger than science fiction. 
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